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The spatial dynamics of pulses in Kerr media with parabolic index profile are examined. It is found that when 
diffraction and gradcd-indcx have opposite signs propagating pulses exhibit an oscillatory pattern, similar to a 
breathing behavior. Furthermore, if the pulse and the index profile are not aligned the pulse oscillates around 
the index origin with frequency that depends on the values of the diffraction and index of refraction. These 
oscillations are not observed when diffraction and gradod-index share the same sign. © 2009 Optical Society 
of America 
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In recent years, great effort has been made in order to ex- 
plain the behavior of light beams propagating through 
interfaces separating optical media with different non- 
linear refractive indices. This interest has been moti- 
vated by different factors that include the analysis of 
pulse propagation and self-phase modulation effects in 
fibers, the close connection with the problem of radi- 
ation mode propagation in the three-dimensional non- 
linear Schrodinger equation, and the persistency in the 
nonlinear wave propagation regime of some properties 
found in the corresponding linear problem, such as the 
occurrence of a parametric instability under modulation 
of parameters [1,2]. 

Nonlinear pulse propagation in gradcd-index optical 
waveguides is a much studied problem, both numeri- 
cally and analytically. In early works, the radial dy- 
namics of pulses in nonlinear fibers were studied by us- 
ing the paraxial ray approximation. However, this was 
shown to give inaccurate results for the self-phase mod- 
ulation of a pulse that propagates in a bulk nonlinear 
medium (see [3] and references therein) and a variational 
method was adopted. Using this approach the authors 
in Ref. [3] studied the resulting dynamics of pulses from 
the combined effects of spatial diffraction, nonlinearity 
and parabolic graded index in radially symmetric fibers. 
More recently, in Ref. [4] the space-time dynamics in 
nonlinear multimodc parabolic index optical fibers were 
studied in the context of a temporal modulation instabil- 
ity induced by spatial transverse effects. We extend this 
approach to include the dynamics of the pulse's center 
in planar waveguides. In fact, we find that the motion 
of the center of the pulse does not depend on any other 
parameter of the pulse and it is only determined by the 
values of the diffraction and the refractive index. 

The nonlinear Schrodinger equation appropriately 
modified to model beam propagation in graded-index, 
nonlinear waveguide amplifiers with refractive index 



no + nix- + 



do 



n2\TpY 



1 



The parameters, do- a-nd n2 are constant, diffraction 
(do) can be positive or negative and the medium can 
be anti-guiding (rii > 0) or guiding (rii < 0). As in 
Ref. [5] , both cases are consider in order to examine the 
interplay between self-defocusing of light and guiding or 
anti-guiding and diffractive effects of the medium. 

Interestingly, Eq. ([T|) is similar to the Gross-Pitaevskii 
equation that describes the dynamics of confined atomic 
Bose-Einstein condensates (BECs) [6,7]. The crucial dif- 
ferences are that, in BECs, do is always positive and 
ni — —fl^ < 0, where is the normalized harmonic 
trap strength. As also shown below, this difference in 
signs is what gives rise to the so-called collective oscilla- 
tions in the BECs context [8]. Furthermore, a dissipative 
variant of this equation has also been used to describe 
the behavior of solitons and self-similar waves in nonlin- 
ear systems exhibiting both spatial inhomogeneity and 
gain or loss at the same time [9, 10]. 

Equation (1) can be restated in variational form using 
the Langrangian density 
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If the injected field is a gaussian beam, in the weakly 
or moderate nonlinear regime {n2 ~ 0) the field remains 
approximately gaussian and the periodic variation of the 
beam parameters along the propagation distance can 
be calculated by means of a variational approach [3,4]. 
Thus, we take 
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where (j){z,x) = a{z)x'^ + (3{z)x + 7(2), describes the 
pulse's phase. For a given value of z, Eq. ([2|) defines a 
gaussian beam invariant along the y direction (planar 
waveguide geometry) displaced by an amount xo along 
the X-axis from the origin of the coordinates. 

Inserting the gaussian ansatz into the expression for 
C and integrating we obtain the average Langrangian of 
the problem, L = J^°° C dx. Using the Euler-Langrange 
equations with variational parameters A, b, xg, a, /? and 
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7 we obtain a set of coupled equations describing the 
evolution of these parameters, namely 



Az = —doaA 
xoz = do(2axo + /3) 
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(3c) 

(3d) 

(3e) 



= \ {b\/2n2A^ - 8dob - A^fln^J^' 



bxn 



+16do&^x5 - 4do/3^) (3f) 

By dividing the first two equations and integrating, we 
obtain A'^ = Eob, where Eq is the energy of the system 
at z = 0. This is equivalent to the conservation of energy 
of the system. Indeed, Eq. ([T]) is integrable and has an 
infinite number of conservation laws the first of which 
characterizes the system's energy, i.e. dx = Eq. 

When we substitute for the profile of Eq. ^ we obtain 
A^ = Eob, as above. In fact, all of Eqs. ^ can also be 
derived using conservation laws arguments. 

In general, the system of Eqs. ([3]) is coupled and non- 
linear. Remarkably, however, differentiating Eq. (|3c[) and 
using the rest of the equations we obtain the uncoupled 
equation 
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The above equation demonstrates that two types of evo- 
lution can be observed. If daiii > the location of the 
center of the pulse moves along an exponential trajec- 
tory. The more interesting case is the one with dofii < 0, 
since it suggests an oscillatory pattern around xq = 0, 
with frequency oj = y/\doni\. This is illustrated in. Fig. 
[U where we show how the center of the pulse evolves, 
with parameters do = —1, ni = ±1.5 (top/bottom) 
and xo(0) = 1. Hereafter, 712 = 1 and the relative 
signs between diffraction and nonlinearity will be con- 
trolled by do. Also, note that depending on the initial 
conditions on xo{z) the oscillation will undergo a sine 
(a;o(0) = 0, Xo,z ^ 0) or cosine (xo(0) ^ 0, xo,^ = 0) 
oscillation. This resembles linear propagation of light in 
a parabolic-index optical fiber, where light is guided if 
the refractive index is maximal in the center, provided 
diffraction is positive. A parabolic index with the min- 
imal value of the index in the center would lead to de- 
focusing of the beam. In the case where light is guided, 
we can expect that different propagation modes will be 
allowed in the waveguiding structure induced by the 
graded index. Therefore, if the input profile is not per- 
fectly symmetric, at least one antisymmetric mode will 
be excited, with a slightly different propagation constant 
than the fundamental mode. The beating between the 
main antisymmetric mode and the main symmetric mode 
will lead to a periodic change in position of the center of 
the beam along the propagation. 
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Fig. 1. (Color Online) The center of a pulse oscillating 
with frequency lo = -^/jdonil, (top) and exponential mov- 
ing (bottom). The contours are the direct evolution of 
Eq. Ill) and the white lines the solution of the uncoupled 
linear equation (width and center of the pulse). 

In a similar manner, a breathing behavior may be ob- 
served if one is to analyze the propagation of a pulsed 
beam through the simple nonlinear device composed of 
a waveguide of certain thickness and refractive index 
n\{^ n\). surrounded by a linear substrate of index n\ 
and a nonlinear cover with Kerr-type nonlinearity of the 
form n = n\ n2|V'P- To demonstrate, we show in Fig. 
[21-(top) the propagation of a unit gaussian under Eq. ([T|) 
with do = ~1 and n\ = 1.5. By changing the values of 
these parameters the qualitative behavior of the oscilla- 
tions are also changing. Moreover, if these two parame- 
ters share the same sign no oscillations are observed, as 
shown in Fig. [21- (bottom). In this case the initial pulse 
(i.e. unit gaussian with equation parameters do ^ —\ 
and n\ = —1.5) is decaying fast. 

To further illustrate this breathing we plot in Fig. [3] 
the evolution of the amplitude A{z) of the pulse for two 
different values of n\ starting from a unit gaussian. The 
analysis for this case is somewhat more involved than in 
the previous case since the equations do not uncouple in 
a trivial way. Indeed, differentiating Eq. (j3dp and using 
the rest of the set we obtain a second order nonlinear, 
coupled equation for q;(z) that reads 
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The oscillatory nature of the equation becomes apparent 
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Fig. 2. (Color Online) The contour plot of the complete 
evolution of a unit gaussian pulse with c?o = — 1, ^t-i = 
1.5 (opposite signs) and a;o = (top). The contour plot 
of the complete evolution of a unit gaussian pulse with 
do = —1, ni = —1.5 (same signs) and xq = (bottom). 
The white line represents the solution of Eqs. dS]). 




Fig. 3. (Color Online) The evolution of the amplitude 
of a unit gaussian with ni ~ —1.5. The diffraction is 
do — ±1 (top/bottom) and and xq = 0. 



from its linear part which can be simplified if 6^ <C 1 to 

O-zz — {idoni)a = 

Again the oscillations exist if doui < and a good esti- 
mate for their frequency is uj' = ^3|c?o'^i I- Notice that 
712 (nonlinearity) comes into the system in higher-order 
and the difference with the frequency of the pulse's cen- 
ter (also apparent in Fig. [4]). 

We finally briefly discuss the evolution of a unit gaus- 
sian originally dislocated at a;o = 1, as shown in Fig. [H 
Since the equation for xo can be uncoupled the oscilla- 
tions of the center of the pulse do not effect the breathing 
of the rest of the pulse's parameters. This means that the 
propagation of the center of the pulse is independent of 
the breathing that may occur due to changes in the re- 
fractive index of the medium. 

To conclude, we demonstrated that under certain con- 
ditions, namely that diffraction and refractive index have 
opposite signs, pulses propagating in parabolic index op- 
tical waveguides exhibit oscillatory patterns. While the 
pattern for the dislocation of the pulse is simple and 
described by a linear equation the rest of the pulse's pa- 
rameters are described by coupled, nonlinear equations. 
In the first case, we provided the exact frequency of the 
oscillation while an estimate based on the linear part of 
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Fig. 4. (Color Online) The complete evolution of a unit 
gaussian centered initially at a;o = 1 and ni = 1.5, do = 
— 1. The white line is the solution of Eqs. ([3]). 
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the equation was provided for the latter. 
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